Meson Supercurrent State in High Density QCD 
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The ground state of three flavor quark matter at asymptotically large density is believed to be the 
color-flavor-locked (CFL) phase. At non-asymptotic density the effect of the non-zero strange quark 
mass cannot be neglected. If the strange quark mass exceeds nis ~ 771^'^ A^^^ the CFL state becomes 
unstable toward the formation of a neutral kaon condensate. Recently, several authors discovered 
that for TTis ~ (2Apf)^''^ the CFL state contains gapless fermions, and that the gapless modes lead 
to an instability in current-current correlation functions. Using an effective theory of the CFL state 
we demonstrate that this instability can be resolved by the formation of a meson supercurrent, 
analogous to Migdal's p-wave pion condensate. This state has a non-zero meson current which is 
canceled by a backflow of gapless fermions. 



Calculations based on weak coupling QCD show that 
the ground state of three flavor QCD at very high baryon 
density is the color-flavor-locked (CFL) phase Q,!^!^. 
During the past several years a lot of effort has been 
devoted to the question how the CFL state is modified 
in the presence of a non-zero strange quark mass. Using 
a BCS-type gap equation it is easy to see that the most 
important parameter is rri^/ (ppA), where is the Fermi 
momentum and A is the gap in the chiral limit 0. 

Using effective field theory methods we showed that 
the CFL state is not modified if rus < m™* ml/'^A'^^^. 
At TOs(crit) the CFL phase undergoes a phase transition 
to a kaon condensed phase Q. This transition is shifted 
to larger values of TOs if instanton effects are important, 
as is the case at moderate density For ml/{2pp) ~ 
2A the effective field theory description breaks down and 
we expect a phase transition to a less symmetric phase 
to occur. More recently, several groups observed that 
gapless fermion modes can appear in the spectrum of the 
CFL phase in the vicinity of the point m'^/{2pp) ~ A 
0, 111 ■ Gapless modes also appear in the kaon condensed 
CFL phase, but the critical strange quark mass is shifted 
to somewhat larger values 

The problem is that gapless fermion modes in a weakly 
coupled (BCS) superfluid tend to cause instabilities in 
current-current correlation functions. These instabilities 
are quite generic and appear in a wide range of sys- 
tems, including the 2 SC p hase, the CFL phase, and cold 
atomic gases 0, Il2j |. The presence of an instabil- 
ity implies that the homogeneous superfluid is not the 
correct ground state 0, Q, ^| . One possibility is that 
the true ground state is an inhomogeneous superconduc- 
tor of the type first considered by Larkin, Ovchinnikov, 
Fulde and Ferrell and generalized to QCD in [Tsllr^. 
In this work we wish to study the possibility that the 
ground state has a non-zero meson current similar to the 
p-wave meson condensate suggested by Migdal, Sawyer 
and Scalapino 0. For this purpose we shall study the 
stability of the s-wave kaon condensate with respect to 
the formation of a non-zero current. Our study is anal- 
ogous to the calculations performed in chiral models of 
nuclear matter |3| and in effective theories of cold atoms 



Our starting point is the effective theory of the CFL 
phase derived in 0, lU 

C = ^Tr (VoSVqS^ - -y^ VSVS) 

+ as (^(Tr [AfS])^ - Tr [AfEAfS] + h.c}j 

+ Tr {NUv^Df^N) - D Tr {nK^'j5 {A^,N}) 
-FTr(iVti;^75[^^,7V]) 

+ ^{{Tr{NN)^[Tr{N)f)+h.c.}. (1) 

The effective lagrangian contains Goldstone boson fields 
E and baryon fields N. The meson fields arise from chiral 
symmetry breaking in the CFL phase T] and the baryon 
fields originate from quark-hadron complementarity |22| . 
The chiral field is given by E = exp(i0''A''//^) where 
/tt is the pion decay constant. M is the quark mass 
matrix. The chiral field and the mass matrix transform 
as E — > LT,R^ and M — > LMR^ under chiral transfor- 
mations (L,i?) e SU{3)l X SU{3)ii. The baryon field 
N transforms according to the adjoint representation of 
flavor SU{3). We will expand N in terms of the baryon 
nonet (p, n, E^, S^, A*'"), = (1, ■[?) is the Fermi ve- 
locity, and A is the superfluid gap. We have suppressed 
the singlet fields associated with the breaking of the exact 
U{l)v and approximate U{1)a symmetries. 

The covariant derivative of the nucleon field is given 
by D^N = d/^iN + «[V^, N]. The vector and axial- vector 
currents are 

v.^-^i^d^^^ + ^^d,^}, A, = -'-^d,j:^) ^, (2) 

where ^ is defined by = E. It follows that ^ transforms 
as ^ ^ L^UixY = U{x)CR^ with U{x) G SU{3)v. 

Symmetry arguments can be used to determine the 
leading mass terms in the effective lagrangian. Be- 
daque and Schafer observed that Xl = MM'^ /{2pp) 
and Xr = M^M/{2pp) enter the microscopic theory like 
the temporal components of left and right handed fla- 
vor gauge fields. We can make the effective lagrangian 
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invariant under this symmetry by introducing the covari- 
ant derivatives 



DqN = doN 



iXr 



(3) 
(4) 



with To = -§{e(ao + iXR)i'^ + eU^o + ^XL)i}■ If the 
density is very large and the QCD couphng is weak then 
the low energy constants can be calculated in perturba- 
tion theory. The leading terms in the meson and baryon 
sector were calculated in Ha and The results are 



f 



21 - 81og(2) 
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as = 



3A2 
47r2 



(5) 



and D ~ F ~ \/2. The gap A can also be computed 
in perturbative QCD, but we will not make use of these 
results and treat A as a parameter. 

The parameter that controls the effect of the strange 
quark mass is the effective chemical potential /is — 
m1/{2pp). When exceeds the mass of the kaon the 
CFL ground state becomes unstable and a Bose conden- 
sate of kaons is formed. If isospin was an exact symmetry 
the energy the if + and meson would be degenerate. 
Because of explicit isospin breaking, and because of the 
constraint of electric charge neutrality, a condensate of 
neutral kaons is favored. The condensed phase is 
characterized by 



1 

£,K0 — I cos(q;/2) isin(a/2) 
isin(a/2) cos(a/2) 



(6) 



with cos(q;) — mj^/^^. From equ. (O we get m| 
3mu{mci + rris) A2/(7r2/2). For typical values of the pa- 
rameters cos(a) ~ (m„/ms)(A//^)2 <^ 1. In the follow- 
ing we shall assume that cos(a) = and the ground state 
is a maximal kaon condensate. 

Properties of the kaon condensed state were studied in 
We showed, in particular, that the kaon condensed 
CFL phase has a gapless fermion mode that appears at 
fig — 4A/3. For comparison, in the ordinary CFL phase 
gapless modes appear at fig = A 0. Recently, several 
groups have shown that gapless fermion modes lead to 
instabilities in the current-current correlation function 
[llL Il2j . Motivated by these results we examine the sta- 
bility of the kaon condensed phase against the formation 
of a non-zero current. 

Consider a spatially varying U{1)y rotation of the 
maximal kaon condensate 
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This state is characterized by non-zero currents 
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In the following we compute the vacuum energy as a func- 
tion of the kaon current Jk = ■ The meson part of 
the effective lagrangian gives a positive contribution 



1 



(10) 



A negative contribution can arise from gapless fermions. 
In order to determine this contribution we have to calcu- 
late the fermion spectrum in the presence of a non-zero 
current. The relevant part of the effective lagrangian is 

£ = Tr {NUv'^Df.N) + Tr (iV^75 (^pA + v ■ Aj 7v) 
A 

-1- — {Tr (NN) - Tr {N) Tr (N) + /i.e.} , (11) 

where we have used D — F = 1/2. The covariant deriva- 
tive is DqN = doN+i[pv, N] and D,N = d,N+iv- [V, TV] 
with V, A given in equ. (|H|) and 



PV,A 



(12) 



The vector potential py and the vector current V are 
diagonal in flavor space while the axial potential pA and 
the axial current A lead to mixing. In order to exhibit 
the mechanism for the instability we will first discuss the 
simple case of vanishing axial vector couplings F = D = 
0. We will then study the realistic case D = F = 1/2. 

We can determine the spectrum by decomposing pv 
and V into isospin and hypcrchargc components. We 
find 

pv = - y/s - -jY, V = — -I3 + —Y, (13) 

where I3 — A3/2 and Y — As/VS are the isospin and hy- 
percharge generators. The dispersion relations can now 
be obtained from the {I3 , Y) quantum numbers of the 
baryon nonet. The result is 



(14) 



v-Jk)- 



where / — v -p — pp is the residual momentum. The neu- 
tral modes (Sqj Aq, Ag) are not affected by the current or 
the vector potential. Equ. (|14|l shows that there are gap- 
less fermion modes as ps approaches 2A, and that gapless 
modes occur earlier if there is a non-zero current. The 
contribution to the vacuum energy from gapless modes 
is 



(15) 
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FIG. 1: Energy density as a function of the current jk for 
several different values of jUs = m'i/{2pF) close to the phase 
transition. 



where the factor 4 is a degeneracy factor and dU, is an 
integral over the Fermi surface. Near the Fermi surface 
we can approximate 



2A 
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V- jk) 



(16) 



The integral in equ. I|15|l receives contributions from one 
of the pole caps on the Fermi surface. The result has 
exactly the same structure as the energy functional of 
a non-relativistic two-component Fermi liquid with non- 
zero polarization which was recently analyzed by Son and 
Stephanov [23|. Introducing dimensionless variables 



Jk , ^J^s- 2A 
aA aA 



(17) 



we can write £ — cAf fh{x) with 



hix) 



X 



(h + xfl'^Qih + x)-{h- xf/^e{h - x) 

(18) 



We have defined the constants 
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M - 



(19) 



where c-^ = (21 — 81og(2))/36 is the numerical coefficient 
that appears in the weak coupling result for /^r. Accord- 
ing to the analysis in [23| the function Jh{x) develops a 
non-trivial minimum if /ii < /i < /12 with hi ~ —0.067 
and /i2 — 0.502. In perturbation theory we find a — 13.9 
and the kaon condensed ground state becomes unstable 
for (A - nJ2) < 0.46A. 

We will see that the instability is weaker if the axial 
coupling is taken into account. The axial vector potential 
and the axial vector current can be written as 



PA 



-v-A = ^ifis + v-jK){e'^''u'^ 



(20) 



where — (Ag ± iXj)/2 are the w-spin raising and 
lowering operators. The axial vector coupling leads 
to mixing among the charged (p, 2~,S^) and neutral 



(n, 5°, S°, A"'®) baryons. Gapless fermions occur in the 
charged sector. The mixing in the charged sector is 



C = 



(21) 

The dispersion relation of the lowest mode is given by 
LOi^A+ ^ (22) 

where lo = (/is + v ■ Jx)/4. The calculation of the energy 
density is completely analogous to the one we performed 
in the previous section. The only difference is that the 
number of gapless modes is reduced, and that the effect 
of the current on the energy of the fermion is smaller. 
The energy functional is £ = aV ft {x) with x defined 
in eq. (|19|l and h = {3 fig — 4A)/(aA). The numerical 
coefficients are c = 2/(15'^c^w^) amd a — 2/(15^cJw^). 

Using the weak coupling values for and Ctt we find 
a — 0.43 and the kaon condensed ground state becomes 
unstable for (A-3/is/4) < 0.007A. In Fig.jTlwe show the 
energy density as a function of the current jk for several 
values of fis near the transition. We note that the mech- 
anism for the instability is quite robust, but the value of 
the onset chemical potential depends sensitively on the 
values of the low energy constants. We also note that 
the homogeneous gapless kaon condensed phase becomes 
stable at even larger values of the effective chemical po- 
tential, (A — 3/Xs/4) < —0.05. This result may not be 
reliable, however, because the contribution from higher 
fermion modes has to be included. 

We note that the ground state has no net current. 
This is clear from the fact that the ground state satisfies 
5£ / (j)K) = 0. As a consequence the meson current 
is canceled by an equal but opposite contribution from 
gapless fermions. This is analogous to what happens in a 
p-wave pion condensate in nuclear matter [T8Lll9l |2^ or 
in the LOFF phase jT^] . We also expect that the ground 
state has no chromo-magnetic instabilities. The effective 
lagrangian given in equ. ^ is formulated in terms of 
gauge invariant fields. Quark-hadron continuity implies 
that we should study the screening length of a SU{3)f 
gauge field. For diagonal gauge fields we find 



0{h) 



(23) 



which shows the magnetic instability for h > and has 
the characteristic square root singularity observed in mi- 
croscopic calculations Tl^. In this work we have calcu- 
lated the energy functional for arbitrary values of jk , not 
just the second derivative at the origin. We find that the 
phase transition to a supercurrent state is first order and 
that the instability sets in for h > he with he < 0- We 
also observe that the second derivative at the new min- 
imum is positive and of the same order of magnitude as 
the screening mass at h = 0. 

Summary: We have that shown that for iig ~ 4A/3, 
which is the threshold for the appearance of gapless 
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fermion modes, the CFL phase becomes unstable toward 
the formation of a meson supercurrent. This state is 
analogous to the p-wave pion condensate studied in low 
density nuclear matter. In that context the name p-wave 
was motivated by the fact that the interaction that leads 
to the formation of a condensate is the v ■ D interaction 
between nucleons and pions that also determines the p- 
wave ttN scattering amplitude. This is also true in our 
case. Supercurrent formation is driven by a v ■ D inter- 
action between baryons and mesons. 

Our result suggests that the series of phases encoun- 
tered in three flavor QCD as the density is lowered start- 
ing from an asymptotically large value consists of the 
CFL phase, an s-wave kaon condensate, and a p-wave 
kaon condensate. There are many questions that remain 
to be addressed. In this work we established the pres- 



ence of an instability but did not determine the exact 
nature of the new ground state. We have to consider 
additional currents, the role of other light fermions, and 
the effects of electric charge neutrality. Once the currents 
become large the effective field theory description breaks 
down and more microscopic approaches along the lines 
of [25I l2fil | are necessary. These methods might also be 
useful in order to determine the fate of the meson super- 
current state as the density is lowered even further and Hs 
approaches 2A. One possibility is that the supercurrent 
phase evolves into a LOFF-type state. 
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